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We present a measurement of Rb, the ratio of the branching fraction for the rare decay D° — > 
K + ir~ to that for the Cabibbo-favored decay D° — > K~ty + . Charge conjugate decays are implicitly 
included. A signal of 2005 ± 104 events for the decay D° —> K + -k~ is obtained using the CDF II 
detector at the Fermilab Tevatron collider. The data set corresponds to an integrated luminosity 
of 0.35 fb _1 produced in pp collisions at ■s/s = 1.96 TeV. Assuming no mixing, we find Rb = 
[4.05 ± 0.21(stat) ±0.11(syst)] x 10 -3 . This measurement is consistent with the world average, and 
comparable in accuracy with the best measurements from other experiments. 

PACS numbers: 13.25.Ft, 14.40.Lb 



The D° can decay to K + ir~ either through a dou- 
bly Cabibbo-suppressed (DCS) tree process or by oscil- 
lation (mixing) to a D° followed by a Cabibbo-favored 
(CF) tree process. The charge-conjugate decays, such as 
£)° — » K~n + , are implied throughout this paper. The 
time-dependent decay rate r(t) for D° — > K + n~ can be 
written in a compact form [l| taking into account the 
experimentally established facts that the rate for mixing 
is at least as small as that for the DCS decay, and that 
the effect of CP violation is small. In this formalism, and 



assuming CP conservation 
r(t) oc e~ rt 



R D + \/iW(rt) + \ v (rt) 2 



(i) 



The parameter Rd is the squared modulus of the ratio 
of DCS to CF amplitudes. The parameters x' and y' 
arc defined in terms of the parameters x — Am/T and 
y = Ar/2r, where Am is the difference in mass between 
the two mass eigenstates, AF is the difference in decay 
width between the two mass eigenstates, and T is the 
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average decay width. The definitions are 

x' = x cos 5 + y sin 6, and y' — —x sin 6 + y cos 5, (2) 

where <5 is the strong phase difference between the DCS 
and CF amplitudes. The mixing rate is small if x 1 -C 1 
and i/'cl. The ratio of branching fractions, 

R B = B(D° ^ K + Tt-)/B(D a ^ K-7T+) (3) 

is given by the ratio of the time-integrals of the corre- 
sponding decay rates, 

x' 2 + v' 2 

R B = R D + ^R^y' + (4) 

Thus, if x' and y' are negligible, the experimentally mea- 
surable quantity R& can be interpreted as the theoretical 
parameter Rd- 

In the limit of flavor SU(3) symmetry, Rd = tan 4 9c 
0,0 where C is the Cabibbo angle, which is measured 
from kaon decays. The world average values Rd = 

(3.62±0.29)x 10~ 3 and tan 4 6 C = (2.88 ±0.27) x 10~ 3 , 
are equal within their uncertainties, consistent with fla- 
vor SU(3) symmetry. However, symmetry violation of 
magnitude less than the current measurement accuracy 
is possible Q if there are differences in the weak decay 
form factors and Fq K , or due to strong interaction 
resonant intermediate states following the charm quark 
decay. 

There is no experimental evidence either for CP viola- 
tion in D° decays or for D°-D° mixing; all measurements 
are consistent with x' = y' = 0. An upper limit on the 
contribution to i?s from mixing can be derived from mea- 
sured upper limits for i',i/'~3x 10~ 2 Q and the world 
average measurement of Rd- A simple estimate, by sub- 
stituting these values into Eq. 10}, gives a contribution 
to Rb less than 2.7 x 1CP 3 , a limit which is comparable 
to the value of Ra- 
in the standard model, theoretical predictions due to 
short-distance weak processes are x',y' ~ 6 x 10 -7 
However, strong interaction effects could result in larger 
values, of order sin 2 9c (=0.048) times an unknown fac- 
tor which describes the size of flavor SU(3) symmetry 
violation. Thus, accurate measurement of Rd and 9c 
can establish the size of the symmetry violation factor 
and make possible the prediction of the standard model 
contribution to mixing. If the standard model contribu- 
tion is small, then D°-D° mixing measurements will be 
sensitive to new physics. Theories involving weak-scale 
supersymmetry or new strong dynamics at the TeV scale 
can accommodate large values of x' and y', up to the 
experimental limits |7J, leaving open the possibility for 
indirect observation of new physics. 

Until now, the most precise measurements of Rd were 
from the B factories. The BABAR collaboration jg re- 
ported R D = [3.59±0.20(stat)±0.27(syst)]x 10~ 3 with a 
DCS signal of 430 events, and the Belle collaboration Q 
reported R D = [3.81 ± 0.17(stat) io;?e( s y st )] x 10 ~ 3 with 
a DCS signal of 845 events. In this paper, using a DCS 



signal of 2005 events and assuming no mixing, we report 
a measurement of Rd with comparable precision to those 
of BABAR and Belle. As in those experiments, we re- 
construct the decay chain D* + — * n + D°, D° — > K + n~, 
where the charge of the 7r + from D* + decay distinguishes 
the D° from its anti-particle, D Q . 

Our measurement uses data collected by the CDF II 
detector at the Fermilab Tevatron collider, from Octo- 
ber 2002 to August 2004. The data corresponds to an 
integrated luminosity of 0.35 fb _1 produced in pp col- 
lisions at ^/s — 1.96 TeV. CDF II is a multi-purpose 
detector with a magnetic spectrometer surrounded by a 
calorimeter and a muon detector. The CDF II compo- 
nents pertinent to this analysis are described briefly be- 
low. A more detailed description is found in ^(j and 
references therein. A cylindrical silicon microstrip vertex 
detector (SVX II) (ll| and a cylindrical drift chamber 
(COT) [la| . immersed in a 1.4 T axial magnetic field, al- 
low reconstruction of tracks (trajectories of charged par- 
ticles) in the pseudo-rapidity range |?7| < 1.3, where 
rj = tanh~ 1 (cos0) and 9 is the angle measured from the 
beamline. The ionization signals from the COT provide 
a measurement of the specific energy loss for a charged 
particle, which is used for particle identification. 

Events were selected in real time using a three-level 
trigger system with requirements developed for a broad 
class of heavy flavor decays. At level 1, tracks are re- 
constructed in the COT in the plane transverse to the 
beamline by a hardware processor (XFT) |T3| . Two op- 
positely charged tracks are required, each with transverse 
momentum greater than 2 GeV/c. In addition, the scalar 
sum of the two transverse momenta must be greater than 
5.5 GeV/c. The opening angle A0 between the two tracks 
in the transverse plane must be less than 135°. At level 2, 
the silicon vertex tracker (SVT) attaches SVX II hits 
to each of the two XFT tracks to increase the measure- 
ment accuracy. The transverse impact parameter do is 
defined as the distance of closest approach, in the trans- 
verse plane, of a track to the beamline. Each of the two 
tracks is required to satisfy 120 ^m < d < 1.0 mm. 
The opening angle cut is tightened (compared to level 1) 
to 2° < Atfi < 90°. The track pair forms a long-lived par- 
ticle candidate which is required to have a decay length 
L xy > 200 /im, where L xy is the transverse distance from 
the beam line to the candidate's vertex, projected along 
the total transverse momentum of the candidate. At level 
3, a conventional computer processor confirms the selec- 
tion with a full event reconstruction. 

The analysis method for determining the ratio of 
branching fractions requires reconstruction of the decay 
chains D*+ -> n+D , D° -> R-TT+ (CF), and D*+ -> 
7r+L>°, D° -> K+n- (DCS). The D° candidate recon- 
struction starts with a pair of oppositely charged tracks 
that satisfy the trigger requirements. The tracks are con- 
sidered with both K~tt + and ir~ K + interpretations. A 
third track, which is required to have px > 0.3 GeV/c, 
is used to form a D* candidate when combined as a 
pion with the D° candidate. The charge of this "tag- 
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ging pion" determines whether the D° candidate decay 
is CF or DCS. 

To reduce systematic uncertainty, the same set of cuts 
is employed for both the CF and DCS decay modes. The 
offline analysis cuts were chosen to maximize the signifi- 
cance of the DCS signal determined from a study of the 
CF signal and the DCS background. The optimization 
was performed without using DCS candidates and before 
the candidates were revealed. The DCS signal was esti- 
mated by scaling the CF signal by the world average for 
Rtj. The DCS background was estimated from candi- 
dates in a control region of D* invariant mass, outside a 
region containing the signal. In the optimization study, 
the same algorithms for data analysis were followed as 
for the DCS signal determination. 

We apply two cuts to reduce the background to the 
DCS signal from CF decays where the D° decay tracks 
are mis-identified. Mis-identification occurs when the 
kaon and pion assignments are mistakenly interchanged. 
This background is characterized by a Ktt mass distribu- 
tion with width about ten times that of the signal peak. 
A DCS candidate that is consistent with being a CF de- 
cay, with K~tt + invariant mass within ± 20 MeV/c 2 of 
the D° mass, is excluded from the DCS signal. This 
cut rejects 97.5% of mis-identified decays, while retain- 
ing 78% of the signal. Since the analysis procedure is the 
same for DCS and CF decays, a CF candidate that is 
consistent with being a DCS decay is excluded from the 
CF signal. 

A cut based on particle identification from specific ion- 
ization in the COT also helps to reject mis-identified de- 
cays, but with a smaller improvement to DCS signal sig- 
nificance than from the cut based on invariant mass. A 
variable Z is defined as the ratio of logarithms of mea- 
sured to predicted charge deposition for a single track. 
The prediction is based on the ionization expected for 
a particle with the measured momentum and a specific 
hypothesis for mass. For a pair of particles, we define 




where the subscripts K and tt indicate the particle hy- 
potheses for the first and second track of the pair, and 
ok-, g-k are the corresponding Gaussian resolutions on 
Z . The particle identification for the pair is chosen by 
the smaller of Sk-k and S^k- This selection is correct 
80.2% of the time, as measured using the CF signal that 
survives the invariant mass cut. 

We apply four cuts to reduce combinatoric background 
from prompt particle production or from improper com- 
binations of tracks in events containing heavy flavor par- 
ticles. These cuts retain most of the signal, with a small 
improvement in the signal significance. Since the D° has 
a long enough mean lifetime to have an observable decay 
length, the decay vertex should be displaced, on aver- 
age, from the production point. We require the trans- 
verse decay length significance L xy /a xy > 5, where 
L xy was defined earlier and a xy is the uncertainty on 
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FIG. 1: Kn invariant mass distribution for candidates re- 
constructed as D° -> K + n~ (DCS), requiring 5 MeV/c 2 < 
Am < 7 MeV/c 2 . The shaded regions are projections from 
the overall fit onto this distribution. This mass plot illustrates 
the relative contributions from the DCS signal and the three 
types of background, as described in the text. 



L xy . The D* has a short enough mean lifetime so that it 
should appear to decay at its production point. The tag- 
ging pion and the D° candidate must be consistent with 
coming from a common point based on a x 2 measure in 
the transverse plane. For the D* vertex, we also require 
\L X y/a xy \ < 15. Furthermore, the tagging pion track 
must have an impact parameter do < 800 /im. 

The Ktt mass distribution for CF decays has been re- 
ported in a recent CDF 11 publication on singly Cabibbo- 
suppressed decays [l5j . The Ktt mass distribution for 
DCS candidates is illustrated by the histogram in Fig. 
for candidates satisfying 5 MeV/c 2 < Am < 7 MeV/c 2 , 
where Am = m{K + TT~TT + ) — m(K + n~) — m(7r + ). Four 
categories of Ktttt combinations contribute to the distri- 
bution. The first category (signal) is DCS signal from 
D*, with the correct D° — > K + tt~ interpretation. The 
second category (random pion) is background from CF 
D° decays, where a randomly selected particle, usually 
from the primary interaction, is used as the tagging pion 
to form the D* candidate. The third category (mis-id 
D°) is background from D* decays where the K and 7r 
assignments from the CF D° decay are mistakenly inter- 
changed. The last category is combinatoric background, 
where one or both tracks do not belong to a D — > Ktt 
decay. 

To determine the signal and background, candidates 
are divided into 60 slices of Am, each slice of width 
0.5 MeV/c 2 . (The distribution in Fig.[T|is for a 2 MeV/c 2 
wide slice for purpose of illustration.) For each slice, 
the Ktt mass distribution is fit using a binned likelihood 
method with a predetermined D° signal shape and a lin- 
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FIG. 2: The number of D° -> K + tv~ (DCS) decays as a 
function of Am. The data points and statistical uncertainty 
bars are taken from the Ktt slice fits. The shaded regions are 
determined from a least-squares fit and show the contributions 
from signal (dark gray) and random tagging pion background 
(light gray) as explained in the text. 



ear function for the background. The D signal shape is 
determined from a fit to the CF Kir distribution, which 
has a negligible background compared to the signal. The 
fit range (1.80 - 1.93) GeV/c 2 excludes background from 
D° — > K + K~ and D° — > Tr + ir~ decays that are mis- 
reconstructed as K + ir~ . The amplitude of the combined 
signal and random pion term and the combinatoric back- 
ground parameters are allowed to float independently for 
each slice. The amplitude and shape of the mis-id D° 
contribution is determined from the CF signal, by inter- 
changing the pion and kaon assignments. 

To determine the amount of DCS signal and random 
pion background, the D° signal yields for the slices are 
plotted as a function of Am, as shown in Fig. [3 This 
distribution is fit using a least-squares method with a sig- 
nal shape predetermined from the CF Am distribution 
and a background function of the form A(Am) B e -c( Am ). 
The amplitudes of the signal and background terms and 
the background shape parameters B and C are deter- 
mined from the fit. The fit results are 2005 ± 104 DCS 
signal and 495172 ± 907 CF signal; their ratio gives 
R B = (4.05 ±0.21) x 10~ 3 . 

Most of the detector properties that affect the DCS 
and CF signals are common and hence do not affect the 
ratio. Thus, there are no systematic uncertainties due to 
geometric acceptance, particle identification, and trigger 
efficiency. While the number of background events is sim- 
ilar for the DCS and CF candidates, the size of the DCS 
signal is much smaller. Thus, systematic uncertainty in 
the DCS background which affects the DCS signal esti- 
mate also affects the ratio. There are three such signifi- 
cant sources of systematic uncertainty, as summarized in 



Table H| 

To estimate the uncertainty due to the assumed com- 
binatoric background shape in the DCS Kir slice fits, 
we compared i?g results for two shapes. The nominal 
shape is linear and gives a good fit. We also tried a 
quadratic form and assigned the change in i?g as a sys- 
tematic uncertainty. To estimate the uncertainty due to 
the assumed Am background shape, we compared Rb 
results for two shapes. The nominal shape is given by 
the function described earlier and gives a good fit. We 
also tried a function with an additional parameter and 
assigned the change in i?g as a systematic uncertainty. 
In fitting the DCS Kir slice fits, the amplitude of the 
mis-id D° background is fixed from the CF signal. A 
simulation of the fitting procedure is used to propagate 
the statistical uncertainty on the background amplitude 
to a systematic uncertainty on R&. 

We considered other sources of systematic uncertainty 
that we found to be negligible. These include effects due 
to small differences in detection efficiencies for K + ver- 
sus K~ and ir + versus ir~ , which are reported in [l5j . 
We tried alternative fits to the DCS Ktt distributions by 
extending the upper limit of the mass range from 1.80 
to 2.00 GeV/c 2 . This study required adding an explicit 
term for background from D° — ► ir + ir~ decays. 

In conclusion, we find R& — [4.05 ± 0.21(stat) ± 
O.ll(syst)] x 10~ 3 . The difference between this value 
and the world average value for tan 4 ^ is (1.17 ± 
0.34) x 10 -3 , a 3.4 a deviation from zero. If not a statis- 
tical fluctuation, this difference could be due to violation 
of flavor SU(3) symmetry causing Rd tan 4 6c, or could 
be a result of mixing. While we cannot rule out the pos- 
sibility of mixing from our result alone, our result is con- 
sistent with the scenario of modest symmetry violation 
and negligible mixing. As shown in Fig.|3 our measured 
value of Rb is in fact consistent with the world average 
and the most accurate individual measurements of Ru 
obtained from BABAR [1] and Belle @- Using the tech- 
nique we have established to extract the D° — > K + ir~ 
signal, we can perform a time-dependent analysis using 
a larger data sample than reported here, to separately 
measure Rd and the mixing parameters x' and y' . 
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TABLE I: Dominant systematic uncertainties for Rb- The 
sources lead to uncertainties in the DCS signal estimate. 
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FIG. 3: Comparison of this measurement of Rd with other 
recent results. All the experimental fits assume no mixing or 
CP violation. The inner set of bars indicate statistical un- 
certainty; the outer set indicates the quadratic sum of statis- 
tical and systematic uncertainties. The shaded region spans 
the PDG average and uncertainty pj. That average includes 
measurements from E791 El, CLEO [13, FOCUS(2001) [H, 
and BABAR 0. The Belle@, FOCUS(2005) [l^ and current 
CDF measurements are not included in the PDG average. 
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